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Abstract 

We propose an entropic geometrical model of crowd behavior dynamics (with dissipa- 
tive crowd kinematics), using Feynman action-amplitude formalism that operates on three 
synergetic levels: macro, meso and micro. The intent is to explain the dynamics of crowds 
simultaneously and consistently across these three levels, in order to characterize their geomet- 
rical properties particularly with respect to behavior regimes and the state changes between 
them. Its most natural statistical descriptor (order parameter) is crowd entropy 5* that sat- 
isfies the Prigogine's extended second law of thermodynamics, dtS > (for any nonisolated 
multi-component system). Qualitative similarities and superpositions between individual and 
crowd configuration manifolds motivate our claim that goal-directed crowd movement operates 
under entropy conservation, dtS = 0, while naturally chaotic crowd dynamics operates under 
(monotonically) increasing entropy function, cits' > 0. Between these two distinct topological 
phases lies a phase transition with a chaotic inter-phase. Both inertial crowd dynamics and 
its dissipative kinematics represent diffusion processes on the crowd manifold governed by the 
Ricci flow. 

Keywords: Crowd behavior dynamics, action-amplitude formalism, entropic crowd mani- 
fold, crowd turbulence, Ricci flow, topological phase transitions. 

1 Introduction 

Today it is well known that disembodied cognition is a myth, albeit one that has had profound 
influence in Western science since Rene Descartes and others gave it credence during the Scientific 
Revolution. In fact, the mind-body separation had much more to do with explanation of method 
than with explanation of the mind and cognition, yet it is with respect to the latter that its impact 
is most widely felt. We find it to be an unsustainable assumption in the realm of crowd behavior. 
Mental intention is (almost immediately) followed by a physical action, that is, a human or animal 
movement |68) . In animals, this physical action would be jumping, running, flying, swimming, 
biting or grabbing. In humans, it can be talking, walking, driving, or kicking, etc. Mathemati- 
cal description of human/animal movement in terms of the corresponding neuro-musculo-skeletal 
equations of motion, for the purpose of prediction and control, is formulated within the realm of 
biodynamics (seeI351|171[3S[371[3Hl[Mllini[lI])- 

The crowd (or, collective) behavior dynamics is clearly formed by some kind of superposition, 
contagion, emergence, or convergence from the individual agents' behavior. According to the 
emergence theory i72 , crowds begin as collectivities composed of people with mixed interests and 
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motives; especially in the case of less stable crowds (expressive, acting and protest crowds) norms 
may be vague and changing; people in crowds make their own rules as they go along. According to 
currently popular convergence theory (see [281 159] ). crowd behavior is not a product of the crowd 
itself, but is carried into the crowd by particular individuals, thus crowds amount to a convergence 
of likeminded individuals. 

We propose that the contagion and convergence theories may be unified by acknowledging that 
both factors may coexist, even within a single scenario: we propose to refer to this third approach 
as behavioral composition. It represents a substantial shift from traditional analytical approaches, 
which have assumed either reduction of a whole into parts or the emergence of the whole from the 
parts. In particular, both contagion and convergence are related to social entropy, which is the 
natural decay of structure (such as law, organization, and convention) in a social system [5]. 

In this paper we attempt to formulate a geometrically predictive model-theory of crowd be- 
havior dynamics, based on the previously formulated individual Life Space Foam concept j46| (see 
Appendix for the brief summary). 

It is today well known that massive crowd movements can be precisely observed/monitored 
from satellites and all that one can see is crowd physics. Therefore, all involved psychology of 
individual crowd agents: cognitive, motivational and emotional - is only a non-transparent in- 
put (a hidden initial switch) for the fully observable crowd physics. In this paper we will label 
this initial switch as 'mental preparation' or 'loading', while the manifested physical action is la- 
beled 'execution'. We propose the entropy formulation of crowd dynamics as a three-step process 
involving individual behavior dynamics and collective behavior dynamics. The chaotic behavior 
phase-transitions embedded in crowd dynamics may give a formal description for a phenomenon 
called crowd turbulence by D. Helbing, depicting crowd disasters caused by the panic stampede 
that can occur at high pedestrian densities and which is a serious concern during mass events like 
soccer championship games or annual pilgrimage in Makkah (see [29 j l30 l [311 154j ). 

2 Generic three— step crowd behavior dynamics 

In this section we propose a generic crowd behavior dynamics as a three-step process based on 
a general partition function formalism. Note that the number of variables Xi in the standard 
partition function from statistical mechanics (see, e.g. [57) ) need not be countable, in which case 
the set of coordinates {x*} becomes a field (j) = 4>{^)- The sum is replaced by the Euclidean path 
integral (that is a Wick-rotated Feynman transition amplitude in imaginary time, see subsection 
[311), as 



More generally, in quantum field theory, instead of the field Hamiltonian H[(j)) we have the action 
S{(t>) of the theory. Both Euclidean path integral. 





real path integral in imaginary time 



(1) 



and Lorentzian one. 




complex path integral in real time 



(2) 
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- represent quantum field theory (QFT) partition functions. We will give formal definitions of the 
above path integrals (i.e., general partition functions) in section 3. For the moment, we only remark 
that the Lorentzian path integral ^ gives a QFT generalization of the (nonlinear) Schrodinger 
equation, while the Euclidean path integral ([T|) in the (rectified) real time represents a statistical 
field theory (SFT) generalization of the Fokker-Planck equation. 

Now, following the framework of the Prigogine's Extended Second Law of Thermodynamics [60j . 
dtS > 0, for entropy S in any complex system described by its partition function, we formulate a 
generic crowd behavior dynamics, based on above partition functions, as the following three-step 
process: 

1. Individual behavior dynamics (I'D) is a transition process from an entropy-growing "load- 
ing" phase of mental preparation, to the entropy-conserving "execution" phase of physical action. 
Formally, XV is given by the phase-transition map: 

"LOADING": atS>0 "EXECUTION": at S=0 



JV : MENTAL PREPARATION =^ PHYSICAL ACTION (3) 
defined by the individual (chaotic) phase-transition amplitude 



/ dtS=o 

/ PHYS. ACTION 



CHAOS 



MENTAL PREP.) := / e^^^i^I*! 
/id J 



where the right-hand-side is the Lorentzian path- integral (or complex path-integral in real time), 
with the individual behavior action 

SiDm = iiD[*]di, 

where -Lid[*I'] is the behavior Lagrangian, consisting of mental cognitive potential and physical 
kinetic energy. 

2. Aggregate behavior dynamics {AT)) represents the behavioral composition-transition map: 

"LOADING" : 9t S>0 "EXECUTION" : dt S=0 



AV : ^ MENTAL PREPARATION =^ ^ PHYSICAL ACTION^ (4) 



ieAD ieAD 



where the (weighted) aggregate sum is taken over all individual agents, assuming equipartition of 
the total behavioral energy. It is defined by the aggregate (chaotic) phase-transition amplitude 



/ dtS=0 

( PHYS. ACTION 



CHAOS 



MENTAl'pREP.) := /p[$]e"^*°[*l 
/ad J 



with the Euclidean path-integral in real time, that is the SFT-partition function, based on the 
aggregate behavior action 

5ad[$]=/ Lj,T,[^]dt, with Lad[*] = 51 
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3. Crowd behavior dynamics (CT)) represents the cumulative transition map: 



"LOADING" : St S>0 "EXECUTION" : dt S=0 



CV : ^ MENTAL PREPARATION => ^ PHYSICAL ACTION, (5) 



where the (weighted) cumulative sum is taken over all individual agents, assuming equipartition 
of the total behavior energy. It is defined by the crowd (chaotic) phase-transition amplitude 



/ dtS=0 

I PHYS. ACTION 



CHAOS 



MENTAL PREP.) :=/P[$]e*^™[ 

/CD J 



with the general Lorcntzian path- integral, that is, the QFT-partition function), based on the crowd 
behavior action 



ScBm^ f Lcorndt, with icD[$] = 51 -^idI*] = -^adI*]- 

ieCD k= # of ADs in CD 



All three entropic phase-transition maps, XT), AD and CD, are spatio-temporal biodynamic 
cognition systems \43 , evolving within their respective configuration manifolds (i.e., sets of their 
respective degrees-of- freedom with equipartition of energy), according to biphasic action-functional 
formalisms with behavior-Lagrangian functions Lid, -^ad and Lqu, each consisting of: 

1. Cognitive mental potential (which is a mental preparation for the physical action), and 

2. Physical kinetic energy (which describes the physical action itself). 

To develop 2D, AD and CD formalisms, we extend into a physical (or, more precisely, biody- 
namic) crowd domain a purely-mental individual Life-Space Foam (LSF) framework for motiva- 
tional cognition [46], based on the quantum-probability concept 



^ The quantum probability concept is based on the following physical facts 1501 1511 

1. The time— dependent Schrodinger equation represents a complex- valued generalization of the real— valued 
Fokker— Planck equation for describing the spatio-temporal probability density function for the system ex- 
hibiting continuous— time Markov stochastic process. 

2. The Feynman path integral (including integration over continuous spectrum and summation over discrete 
spectrum) is a generalization of the time— dependent Schrodinger equation, including both continuous-time 
and discrete-time Markov stochastic processes. 

3. Both Schrodinger equation and path integral give 'physical description' of any system they are modelling in 
terms of its physical energy, instead of an abstract probabilistic description of the Fokker-Planck equation. 

Therefore, the Feynman path integral, as a generalization of the (nonlinear) time— dependent Schrodinger equation, 
gives a unique physical description for the general Markov stochastic process, in terms of the physically based 
generalized probability density functions, valid both for continuous-time and discrete— time Markov systems. Its 
basic consequence is this: a different way for calculating probabilities. The difference is rooted in the fact that sum 
of squares is different from the square of sums, as is explained in the following text. Namely, in Dirac-Feynman 
quantum formalism, each possible route from the initial system state A to the final system state B is called a history. 
This history comprises any kind of a route, ranging from continuous and smooth deterministic (mechanical-like) 
paths to completely discontinues and random Markov chains (see, e.g., 1151 ). Each history (labelled by index i) is 
quantitatively described by a complex number. 

In this way, the overall probability of the system's transition from some initial state A to some final state B is 
given not by adding up the probabilities for each history-route, but by 'head— to-tail' adding up the sequence of 
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The behavioral approach to TD, AT) and CD is based on entropic motor control [331 I34j . 
which deals with neuro-physiological feedback information and environmental uncertainty. The 
probabilistic nature of human motor action can be characterized by entropies at the level of the 
organism, task, and environment. Systematic changes in motor adaptation are characterized as 
task-organism and environment-organism tradeoffs in entropy. Such compensatory adaptations 
lead to a view of goal-directed motor control as the product of an underlying conservation of entropy 
across the task-organism-environment system. In particular, an experiment conducted in [34] 
examined the changes in entropy of the coordination of isometric force output under different levels 
of task demands and feedback from the environment. The goal of the study was to examine the 
hypothesis that human motor adaptation can be characterized as a process of entropy conservation 
that is reflected in the compensation of entropy between the task, organism motor output, and 
environment. Information entropy of the coordination dynamics relative phase of the motor output 
was made conditional on the idealized situation of human movement, for which the goal was always 
achieved. Conditional entropy of the motor output decreased as the error tolerance and feedback 
frequency were decreased. Thus, as the likelihood of meeting the task demands was decreased 
increased task entropy and/or the amount of information from the environment is reduced increased 
environmental entropy, the subjects of this experiment employed fewer coordination patterns in 
the force output to achieve the goal. The conservation of entropy supports the view that context 
dependent adaptations in human goal-directed action are guided fundamentally by natural law 
and provides a novel means of examining human motor behavior. This is fundamentally related to 
the Heisenberg uncertainty principle |5T] and further supports the argument for the primacy of a 
probabilistic approach toward the study of biodynamic cognition systems!^ 

amplitudes making-up each route first (i.e., performing the sum— over— histories) - to get the total amplitude as a 
'resultant vector', and then squaring the total amplitude to get the overall transition probability. 

Here we emphasize that the domain of validity of the 'quantum' is not restricted to the microscopic world 1731 . 
There are macroscopic features of classically behaving systems, which cannot be explained without recourse to the 
quantum dynamics. This field theoretic model leads to the view of the phase transition as a condensation that is 
comparable to the formation of fog and rain drops from water vapor, and that might serve to model both the gamma 
and beta phase transitions. According to such a model, the production of activity with long-range correlation in 
the brain takes place through the mechanism of spontaneous breakdown of symmetry (SBS), which has for decades 
been shown to describe long-range correlation in condensed matter physics. The adoption of such a field theoretic 
approach enables modelling of the whole cerebral hemisphere and its hierarchy of components down to the atomic 
level as a fully integrated macroscopic quantum system, namely as a macroscopic system which is a quantum system 
not in the trivial sense that it is made, like all existing matter, by quantum components such as atoms and molecules, 
but in the sense that some of its macroscopic properties can best be described with recourse to quantum dynamics 
(see 1141 and references therein). Also, according to Freeman and Vitielo, many-body quantum field theory appears 
to be the only existing theoretical tool capable to explain the dynamic origin of long-range correlations, their rapid 
and efficient formation and dissolution, their interim stability in ground states, the multiplicity of coexisting and 
possibly non-interfering ground states, their degree of ordering, and their rich textures relating to sensory and motor 
facets of behaviors. It is historical fact that many-body quantum field theory has been devised and constructed 
in past decades exactly to understand features like ordered pattern formation and phase transitions in condensed 
matter physics that could not be understood in classical physics, similar to those in the brain. 

^Our entropic action-amplitude formalism represents a kind of a generalization of the Haken-Kelso-Bunz (HKB) 
model of self-organization in the individual's motor system |16l I55| . including: multi-stability, phase transitions 
and hysteresis effects, presenting a contrary view to the purely feedback driven systems. HKB uses the concepts 
of synergetics (order parameters, control parameters, instability, etc) and the mathematical tools of nonlinearly 
coupled (nonlinear) dynamical systems to account for self-organized behavior both at the cooperative, coordinative 
level and at the level of the individual coordinating elements. The HKB model stands as a building block upon 
which numerous extensions and elaborations have been constructed. In particular, it has been possible to derive 
it from a realistic model of the cortical sheet in which neural areas undergo a reorganization that is mediated by 
intra- and inter-cortical connections. Also, the HKB model describes phase transitions ('switches') in coordinated 
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Yet, it is well known that humans possess more degrees of freedom than are needed to perform 
any defined motor task, but are required to co-ordinate them in order to reliably accomplish high- 
level goals, while faced with intense motor variability. In an attempt to explain how this takes place, 
Todorov and Jordan have formulated an alternative theory of human motor co-ordination based on 
the concept of stochastic optimal feedback control [TO]- They were able to conciliate the require- 
ment of goal achievement (e.g., grasping an object) with that of motor variability (biomechanical 
degrees of freedom). Moreover, their theory accommodates the idea that the human motor control 
mechanism uses internal 'functional synergies' to regulate task-irrelevant (redundant) movement. 

Also, a developing field in coordination dynamics involves the theory of social coordination, 
which attempts to relate the DC to normal human development of complex social cues following 
certain patterns of interaction. This work is aimed at understanding how human social interaction 
is mediated by meta-stability of neural networks [49]. fMRI and EEG are particularly useful in 
mapping thalamocortical response to social cues in experimental studies. In particular, a new 
theory called the Phi complex has been developed by S. Kelso and collaborators, to provide exper- 
imental results for the theory of social coordination dynamics (see the recent nonlinear dynamics 
paper discussing social coordination and EEG dynamics [7T]). According to this theory, a pair 
of phi rhythms, likely generated in the mirror neuron system, is the hallmark of human social 
coordination. Using a dual-EEG recording system, the authors monitored the interactions of eight 
pairs of subjects as they moved their fingers with and without a view of the other individual in 
the pair. 

3 Formal model of crowd dynamics 

In this section we formally develop a three-step crowd behavior dynamics, conceptualized by tran- 
sition maps (I3])-(|3])-(l5]), in agreement with Haken's synergetics [T71[TH]. We first develop a macro- 
level individual behavior dynamics IP. Then we generalize XV into an 'orchestrated' behavioral- 
compositional crowd dynamics CD, using a quantum-like micro-level formalism with individual 
agents representing 'crowd quanta'. Finally we develop a meso-level aggregate statistical-field 
dynamics AV, such that composition of the aggregates AT> makes-up the crowd. 

3.1 Individual behavior dynamics (I'D) 

I'D transition map ([3]) is developed using the following action-amplitude formalism (see [46|, I45j): 
1. Macroscopically, as a smooth Riemannian n— manifold Afjo with steady force-fields and 

human movement as follows: (i) when the agent begins in the anti-phase mode and speed of movement is increased, 
a spontaneous switch to symmetrical, in-phase movement occurs; (ii) this transition happens swiftly at a certain 
critical frequency; (iii) after the switch has occurred and the movement rate is now decreased the subject remains 
in the symmetrical mode, i.e. she does not switch back; and (iv) no such transitions occur if the subject begins 
with symmetrical, in-phase movements. The HKB dynamics of the order parameter relative phase as is given by a 
nonlinear first-order ODE: 

<j>= {a + 2l3r^) sin<j> - f3r^ sin2(p, 

where </> is the phase relation (that characterizes the observed patterns of behavior, changes abruptly at the transition 
and is only weakly dependent on parameters outside the phase transition), r is the oscillator amplitude, while a, 13 
are coupling parameters (from which the critical frequency where the phase transition occurs can be calculated). 
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behavioral paths, modelled by a real-valued classical action functional 6*10 [^'], of the form 



Sium = J^'"^ Liumdt, 

(where macroscopic paths, fields and geometries are commonly denoted by an abstract field 
symbol $') with the potential-energy based Lagrangian L given by 

where C is Lagrangian density, the integral is taken over all n local coordinates = (t) of 
the ID, and d^j^'' are time and space partial derivatives of the variables over coordinates. 
The standard least action principle 

<55id[$] =0, 

gives, in the form of the Euler-Lagrangian equations, a shortest path, an extreme force-field, 
with a geometry of minimal curvature and topology without holes. We will see below that 
high Riemannian curvature generates chaotic behavior, while holes in the manifold produce 
topologically induced phase transitions. 

2. Microscopically, as a collection of wildly fluctuating and jumping paths (histories), force-fields 
and geometries/topologies, modelled by a complex-valued adaptive path integral, formulated 
by defining a multi-phase and multi-path (multi-field and multi-geometry) transition ampli- 
tude from the entropy-growing state of Mental Preparation to the entropy-conserving state 
of Physical Action, 

(Physical Action I Mental Preparation)iD / e*'^™!*! (6) 

JlD 

where the functional ID-measure I?[w<I'] is defined as a weighted product 

N 

2?[w<i>] = lim TTwgd*^, {i ^ 1, ...,n = con + dis), (7) 

s=l 

representing an 00— dimensional neural network 46], with weights Ws updating by the general 
rule 

new value{t + 1) — old value{t) + innovation{t). 

More precisely, the weights Ws — Ws{t) in ([54|) are updated according to one of the two 
standard neural learning schemes, in which the micro-time level is traversed in discrete 
steps, i.e., li t = t^^ti^ ...^ts then t + 1 = ti, t2, ts+iH 

■^The traditional neural networks approaches are known for their classes of functions they can represent. Here 
we are talking about functions in an extensional rather than merely intensional sense; that is, function can be read 
as input/output behavior ISlBl [mi26l . This limitation has been attributed to their low-dimensionality (the largest 
neural networks are limited to the order of 10^ dimensions [53]). The proposed path integral approach represents a 
new family of function-representation methods, which potentially offers a basis for a fundamentally more expansive 
solution. 
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(a) A self-organized, unsupervised (e.g., Hebbian-like [27j ) learning rule: 

ws{t+l)=ws{t) + ^{wt{t)~w:{t)), (8) 

where a — a{t), rj — r]{t) denote signal and noise, respectively, while superscripts d and 
a denote desired and achieved micro-states, respectively; or 

(b) A certain form of a supervised gradient descent learning: 

w,{t+l) = w,{t) -r]\7J{t), (9) 

where 77 is a small constant, called the step size, or the learning rate, and V J(n) denotes 
the gradient of the 'performance hyper-surface' at the <— th iteration. 

(Note that we could also use a reward-based, reinforcement learning rule [53], in which 
system learns its optimal policy: innovation{t) = \reward(t) — penalty (t)\. ) 

In this way, we effectively derive a unique and globally smooth, causal and entropic phase-transition 
map ^ , performed at a macroscopic (global) time-level from some initial time tini to the final time 
tfin. Thus, we have obtained macro-objects in the ID: a single path described by Newtonian-like 
equation of motion, a single force-field described by Maxwellian-like field equations, and a single 
obstacle-free Riemannian geometry (with global topology without holes) . 

In particular, on the macro-level, we have the ID-paths, that is biodynamical trajectories 
generated by the Hamilton action principle 

(55idN = 0, 

with the Newtonian action 6*10 [a;] given by (Einstein's summation convention over repeated indices 
is always assumed) 

S'idN = J [ip + -gij x'i^]dt, (10) 
where = (p{t, a;*) denotes the mental LSF-potential field, while the second term, 

represents the physical (biodynamic) kinetic energy generated by the Riemannian inertial metric 
tensor gij of the configuration biodynamic manifold Mid (see Figure [1]). The corresponding Euler- 
Lagrangian equations give the Newtonian equations of human movement 

jT^^-T^.^F,, (11) 

where subscripts denote the partial derivatives and we have defined the covariant muscular forces 
Fi = Fi{t, x\ i*) as negative gradients of the mental potential (p{x'^), 

F,^~ip,.. (12) 

Equation (jlip can be put into the standard Lagrangian form as 

j^L^.^L.,., with L = T-^{x'), (13) 
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or (using the Legendre transform) into the forced, dissipative Hamiltonian form [361 139] 

x'^dpM + dp^R, F,-d^.H + d^^R, (14) 

where pi are the generaUzed momenta (canonically-conjugate to the coordinates x*), H = H(p, x) 
is the Hamiltonian (total energy function) and R = R{p, x) is the general dissipative function. 




Figure 1: Riemannian configuration manifold Mid of human biodynamics is defined as a topological 
product M = Yli SE{3y of constrained Euclidean 5i?(3)-groups of rigid body motion in 3D Euclidean 
space (see [411144) 1. acting in all major (synovial) human joints. The manifold Af is a dynamical structure 
activated/controlled by potential covariant forces (|12|l produced by a synergetic action of about 640 skeletal 
muscles [39) . 

The human motor system possesses many independently controllable components that often 
allow for more than a single movement pattern to be performed in order to achieve a goal. Hence, 
the motor system is endowed with a high level of adaptability to different tasks and also environ- 
mental contexts |34) . The multiple SE(3)-dynamics applied to human musculo-skeletal system 
gives the fundamental law of biodynamics, which is the covariant force ZawQ 

Force co-vector field — Mass distribution x Acceleration vector-field, (15) 

which is formally written: 

Ft = grja\ = l,...,n = dim(M)), 

where Fi are the covariant force/torque components, gij is the inertial metric tensor of the con- 
figuration Riemannian manifold M = J^^ SE{3y {gij defines the mass-distribution of the human 
body), while are the contravariant components of the linear and angular acceleration vector- 
field. Both Lagrangian and (topologically equivalent) Hamiltonian development of the covariant 

*This fundamental biodynamic law states that contrary to common perception, acceleration and force are not 
quantities of the same nature: while acceleration is a non-inertial vector-field, force is an inertial co-vector-field. 
This apparently insignificant difference becomes crucial in injury prediction/prevention, especially in its derivative 
form in which the 'massless jerk' (= d) is relatively benign, while the 'massive jolt' (= F) is deadly. 
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force law is fully elaborated in [3^1 SOI SIl HI] • This is consistent with the postulation that human 
action is guided primarily by natural law [56] , 

On the micro-ID level, instead of each single trajectory defined by the Newtonian equation of 
motion (jlip . we have an ensemble of fluctuating and crossing paths on the configuration manifold 
M with weighted probabilities (of the unit total sum). This ensemble of micro-paths is defined by 
the simplest instance of our adaptive path integral ^ , similar to the Feynman's original sum over 
histories, 

(Physical Action I Mental Preparation) Af = / 'D[wx\e'^^'^\ (16) 

where 'D[wx\ is the functional ID-measure on the space of all weighted paths, and the exponential 
depends on the action 6*10 [2;] given by ([TU)) . 



3.2 Crowd behavioral— compositional dynamics (CP) 



In this subsection we develop a generic crowd CD, as a unique and globally smooth, causal and 
entropic phase-transition map ([5]), in which agents (or, crowd's individual entities) can be both 
humans and robots. This crowd behavior action takes place in a crowd smooth Riemannian 
3n— manifold M . Recall from Figure[l]that each individual segment of a human body moves in the 
Euclidean 3-space according to its own constrained SE(3)-group. Similarly, each individual 
agent's trajectory, = x^{t), i = 1, ■■■n, is governed by the Euclidean SE(2)-group of rigid body 
motions in the plane. (Recall that a Lie group SE{2) = SO (2) x R is a set of all 3 x 3— matrices 
of the form: 

cos 6* sin 6* x 



sm f 




cost 




including both rigid translations (i.e., Cartesian x, ?/— coordinates) and rotation matrix 



cos f 
— sin( 



sm t 
cos( 



in Euclidean plane (see [lUHll).) The crowd configuration manifold M is defined as a union 
of Euclidean SE(2)-groups for all n individual agents in the crowd, that is crowd's configuration 
3n— manifold is defined as a set 



M = 



k=l 



SE{2f = ^SO{2f 



(17) 



k=l 



coordinated by x''' = {a;'', y'', 6'''}, (for fc = 1, 2, n). 



In other words, the crowd configuration manifold M is a dynamical planar graph with individual 
agents' SE(2)-groups of motion in the vertices and time-dependent inter-agent distances lij = 
\x'^{ti) — x^{tj)\ as edges. 

Similarly to the individual case, the crowd action functional includes mental cognitive potential 
and physical kinetic energy, formally given by (with i, j = 1, ...,3n): 



^ [x , X'^ ,ti,tj^ — 

with 



1 
2 



5{ltj) x\U)i3{tj) dUdtj + - g^jx'{t)x^t)dt 



[x' {U ) - {tj )] ^ , where IN <U,tj,t< OUT. 



(18) 
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The first term in represents the mental potential for the interaction between any two agents 
and within the total crowd matrix x'^. (Although, formally, this term contains cognitive 
velocities, it still represents 'potential energy' from the physical point of view.) It is defined as a 
double integral over a delta function of the square of interval between two points on the paths 
in their individual cognitive LSFs. Interaction occurs only when this LSF-distance between the 
two agents cc* and vanishes. Note that the cognitive intentions of any two agents generally 
occur at different times ti and tj unless ti = tj, when cognitive synchronization occurs. This term 
effectively represents the crowd cognitive controller (see [45j ) . 

The second term in p8l) represents kinetic energy of the physical interaction of agents. Namely, 
after the above cognitive synchronization is completed, the second term of physical kinetic energy 
is activated in the common CD manifold, reducing it to just one of the agents' individual manifolds, 
which is equivalent to the center-of-mass segment in the human musculo-skeletal system. Therefore, 
from we can derive a generic Euler-Lagrangian dynamics that is a composition of p^ . which 
also means that we have in place a generic Hamiltonian dynamics that is a amalgamate of ()14p . 
and the crowd covariant force law ()15|) . the governing law of crowd biodynamics: 

Crowd force co- vector field = Crowd mass distribution x Crowd acceleration vector-field, 
formally: Fi — gija\ where gij is the inertial metric tensor of crowd manifold Af. (19) 

The left-hand side of this equation defines forces acting on the crowd, while right-hand defines its 
mass distribution coupled to the crowd kinematics (C/C, described in the next subsection). 

At the slave level, the adaptive path integral, representing an oo— dimensional neural network, 
corresponding to the crowd behavior action (fTS)) . reads 

(Physical Action | Mental Preparation) cd = / ^?[«^, x, y] e*'4[^'f'*''*jl, (20) 

where the Lebesgue-type integration is performed over all continuous paths — x'^{ti) and — 
y-'itj), while summation is performed over all associated discrete Markov fluctuations and jumps. 
The symbolic differential in the path integral ((20)) represents an adaptive path measure, defined as 
the weighted product 

N 

T>[w,x,y] = \im Tlw^^.dx'dy^, (i, j = 1, n). (21) 

N^oo •' 
s=l 

The quantum-field path integral (|20|) - pT|) defines the microstate CP— level, an ensemble of fluc- 
tuating and crossing paths on the crowd 3n— manifold M . 

3.3 Dissipative crowd kinematics (C/C) 

The crowd action (jlSp with its amalgamate Lagrangian dynamics (jl3p and amalgamate Hamilto- 
nian dynamics p4p . as well as the crowd force law (jl9p define the macroscopic crowd dynamics, 
CD. Suppose, for a moment, that CD is force-free and dissipation free, therefore conservative. 
Now, the basic characteristic of the conservative Lagrangian/Hamiltonian systems evolving in the 
phase space spanned by the system coordinates and their velocities/momenta, is that their flow 
iff' (explained below) preserves the phase-space volume. This is proposed by the Liouvillc theo- 
rem, which is the well known fact in statistical mechanics. However, the preservation of the phase 



11 



volume causes structural instability of the conservative system, i.e., the phase-space spreading 
effect by which small phase regions Rt will tend to get distorted from the initial one Ro during 
the conservative system evolution. This problem, governed by entropy growth {dtS > 0), is much 
more serious in higher dimensions than in lower dimensions, since there are so many 'directions' 
in which the region can locally spread (see [631 141j ). This phenomenon is related to conservative 
Hamiltonian chaos (see section 2] below) . 

However, this situation is not very frequent in case of 'organized' human crowd. Its self- 
organization mechanisms are clearly much stronger than the conservative statistical mechanics 
effects, which we interpret in terms of Prigogine's dissipative structures. Formally, if dissipa- 
tion of energy in a system is much stronger then its inertial characteristics, then instead of the 
second-order Newton-Lagrangian dynamic equations of motion, we are actually dealing with the 
first-order driftless (non-acceleration, non-inertial) kinematic equations of motion, which is related 
to dissipative chaos [61]. Briefly, the dissipative crowd flow can be depicted like this: from the set 
of initial conditions for individual agents, the crowd evolves in time towards the set of the corre- 
sponding entangled attractors^ which are mutually separated by fractal (non-integer dimension) 
separatrices. 

In this subsection we elaborate on the dissipative crowd kinematics (C/C), which is self-controlled 
and dominates the CD if the crowd's inertial forces are much weaker then the the crowd's dissipation 
of energy, presented here in the form of nonlinear velocity controllers. 

Recall that the essential concept in dynamical systems theory is the notion of a vector-field (that 
we will denote by a boldface symbol), which assigns a tangent vector to each point p in the manifold 
in case. In particular, v is a gradient vector-field if it equals the gradient of some scalar function. 
A flow-line of a vector-field v is a path ^{t) satisfying the vector ODE, j{t) = v{'j{t)), that is, 
V yields the velocity field of the path 7(t). The set of all flow lines of a vector-field v comprises 
its flow (fit that is (technically, see e.g., [HI [44]) a one-parameter Lie group of diffeomorphisms 

^Recall that quantum entanglement is a quantum mechanical phenomenon in which the quantum states of two 
or more objects are linked together so that one object can no longer be adequately described without full mention 
of its counterpart — even though the individual objects may be spatially separated. This interconnection leads to 
correlations between observable physical properties of remote systems. The related phenomenon of wave-function 
collapse gives an impression that measurements performed on one system instantaneously influence the other systems 
entangled with the measured system, even when far apart. 

Entanglement has many applications in quantum information theory. Mixed state entanglement can be viewed 
as a resource for quantum communication. A common measure of entanglement is the entropy of a mixed quantum 
state (see, e.g. |51)). Since a mixed quantum state p is a probability distribution over a quantum ensemble, this leads 
naturally to the definition of the von Neumann entropy, S{p) = — Tr (p log2 p) , which is obviously similar to the 
classical S/iannon eniropi/ for probability distributions (pi ,•■ ■ , pn ), defined as 5(pi,--- , p?i) = — P» l°g2 P»- 
in statistical mechanics, one can say that the more uncertainty (number of microstates) the system should possess, 
the larger is its entropy. Entropy gives a tool which can be used to quantify entanglement. If the overall system is 
pure, the entropy of one subsystem can be used to measure its degree of entanglement with the other subsystems. 

The most popular issue in a research on dissipative quantum brain modelling has been quantum entanglement 
between the brain and its environment |65l 1661 . where the brain— environment system has an entangled 'memory' 
state, identified with the ground (vacuum) state |0 >Ar, that cannot be factorized into two single-mode states. 
(In the Vitiello-Pessa dissipative quantum brain model '65' '66J, the evolution of the A/"- coded memory system was 
represented as a trajectory of given initial condition running over time-dependent states \0{t) >j\f, each one mini- 
mizing the free energy functional.) Similar to this microscopic brain-environment entanglement, we propose a kind 
of macroscopic entanglement between the operating modes of the crowd behavior controller and its biodynamics, 
which can be considered as a 'long-range correlation'. 

Applied externally to the dimension of the crowd 3n— manifold M, entanglement effectively reduces the number 
of active degrees of freedom in I I17I I. 
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(smooth bijective functions) generated by a vector-field v on M, such that 

tpt o tps = ipt+s, (^0 = identity, which gives: ^{t) = ipt{'y{0)). 

Analytically, a vector-field v is defined as a set of autonomous ODEs. Its solution gives the flow 

Lpt, consisting of integral curves (or, flow lines) j{t) of the vector-fleld, such that all the vectors 
from the vector-field are tangent to integral curves at different representative points p € M. In this 
way, through every representative point p € M passes both a curve from the flow and its tangent 
vector from the vector-field. Geometrically, vector-field is defined as a cross-section of the tangent 
bundle TM of the manifold M. 

In general, given an nD frame {dt} = {d/dx^} on a smooth n— manifold M (that is, a basis of 
tangent vectors in a local coordinate chart a;' = {x^, ...,x") C M), we can define any vector-field 
V on M by its components = v^{t) as 

^ _ ^ig _ ^ _ ^1 ^ _|_ 9 
' dx^ dx^ "' dx"' 

Thus, a vector-field v e X{M) (where X{M) is the set of all smooth vector-fields on M) is actually 
a differential operator that can be used to differentiate any smooth scalar fimction / = f{x^, .t") 
on M, as a directional derivative of / in the direction of v. This is denoted simply v/, such that 

In particular, if v = 7(f) is a velocity vector-field of a space curve ^{t) = {x^{t), ...,x"{t)), 
defined by its components = a;'(t), directional derivative of /(a;*) in the direction of v becomes 

which is a rate-of-change of / along the curve ^{t) at a point a;*(t). 

Given two vector-fields, u = u^di,v = v^di G X{M), their Lie bracket (or, commutator) is 
another vector-field [u, v] e X{M), defined by 

[u, v] = uv — vu = v^dj — v^dj u^di, 

which, applied to any smooth function / on M, gives 

[u,v](/)=u(v(/))-v(u(/)). 

The Lie bracket measures the failure of 'mixed directional derivatives' to commute. Clearly, 
mixed partial derivatives do commute, dj\ — 0, while in general it is not the case, [u, v] ^ 0. 
In addition, suppose that u generates the flow ipt and v generates the flow ips- Then, for any 
smooth function / on M, we have at any point p on M, 

[U, v](/)(p) = ^ (/(<Pa(<^t(p))) - fivMp))), 

which means that in f{ips{(ptip)) we are starting at p, flowing along v a little bit, then along u 
a little bit, and then evaluating /, while in /{^ptifsip)) we are flowing first along u and then v. 
Therefore, the Lie bracket infinitesimally measures how these fiows fail to commute. 
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The Lie bracket satisfies the fohowing three properties (for any three vector-fields u, v, w G M 
and two constants a,b - thus forming a Lie algebra on the crowd manifold M): 

(i) [u, v] = — [v, u]— skew-symmetry; 

(ii) [u, av + bw] = a[u, v] + b[u, w]— bilinearity; and 

(iii) [u, [v, w]] -I- [v, [w, u]] -I- [w, [u, v]]— Jacobi identity. 

A new set of vector-fields on M can be generated by repeated Lie brackets of u, v, vif G M. 

The Lie bracket is a standard tool in geometric nonlinear control theory (see, e.g. [SI [44]). Its 
action on vector-fields can be best visualized using the popular car parking example, in which the 
driver has two different vector-field transformations at his disposal. They can turn the steering 
wheel, or they can drive the car forward or backward. Here, we specify the state of a car by 
four coordinates: the {x, y) coordinates of the center of the rear axle, the direction 6 of the 
car, and the angle (f> between the front wheels and the direction of the car. I is the constant 
length of the car. Therefore, the 4D configuration manifold of a car is a set M = 50(2) x 
R^, coordinated by x = {x, y, 9, (j>}, which is slightly more complicated than the individual crowd 
agent's 3D configuration manifold SE{2) = SO{2) x R, coordinated by x = {x, y, 9}. The driftless 
car kinematics can be defined as a vector ODE: 



X = U(x) Cl + V(x) C2, 



(22) 



with two vector-fields, u, v G X(M), and two scalar control inputs, ci and C2. The infinitesimal 
car-parking transformations will be the following vector-fields 



/ ^ n d . „ d 

u(x) = DRIVE = cos b*- — hsmfc'— 

ox oy 



and 



v(x) 



STEER = 



/ \ 




V 1 J 



tan0 d 
I 89 



The car kinematics 



therefore expands into a matrix ODE: 



/ cos 6* 
sin 9 
J tan ( 




/ x\ 

y 

9 



DRIVE • Cl -I- STEER • C2 





/ cos6i \ 








sin6' 







C2 = 


J tan0 


•Cl + 







I ) 




1 1 / 



C2 . 



However, steer and drive do not commute (otherwise we could do all your steering at home 
before driving of on a trip). Their combination is given by the Lie bracket 



[v, Uj = [STEER, DRIVEj = 



1 



d 



I cos^ <j) 89 



= WRIGGLE. 



The operation [v, u] = wriggle = [steer, drive] is the infinitesimal version of the sequence of 
transformations: steer, drive, steer back, and drive back, i.e., 

{steer, drive, steer"^, drive"^}. 
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Now, WRIGGLE can get us out of some parking spaces, but not tight ones: we may not have enough 
room to WRIGGLE out. The usual tight parking space restricts the drive transformation, but not 
STEER. A truly tight parking space restricts steer as well by putting your front wheels against 
the curb. 

Fortunately, there is still another commutator available: 

[u, [v, u]] = [drive, [steer, drive]] = [[u, v], u] = 

DRIVE, wriggle = — SUl9- COS^— = SLIDE. 

/ cos^ (f) \ OX oy J 

The operation [[u, v],u] = slide = [drive, wriggle] is a displacement at right angles to the car, 
and can get us out of any parking place. We just need to remember to steer, drive, steer back, 
drive some more, steer, drive back, steer back, and drive back: 

{steer, drive, steer~\ drive, steer, drive^\ steer^\ drive^^}. 

We have to reverse steer in the middle of the parking place. This is not intuitive, and no doubt is 
part of a common problem with parallel parking. 

Thus, from only two controls, ci and C2, we can form the vector-fields drive = u, steer = v, 
WRIGGLE = [v, u], and slide = [[u, v],u], allowing us to move anywhere in the car configuration 
manifold M = SO(2) x K.^. All above computations are straightforward in Mathematical '^'^ if 
we define the following three symbolic functions: 

1. Jacobian matrix: JacMat[vXist, x_List] := Outer[D, v, x]; 

2. Lie bracket: LieBrc[u_List, v_List, x_List] :— JacMat[v, x] . u - JacMat[u, x] . v; 

3. Repeated Lie bracket: Adj[u_List, v_List, x_List, k_] :— 

If[k —— 0, V, LieBrc[u, Adj[u, v, x, k - 1], x] ]; 
In case of the human crowd, we have a slightly simpler, but multiplied problem, i.e., super- 
position of n individual agents' motions. So, we can define the dissipative crowd kinematics as a 
system of n vector ODEs: 

x*^' = u'=(x) + v'=(x) 4, where (23) 

cos*^ ( 

u''(x) = drive''' = cos*" 6*-^ + sin*^ = I sin*" 6* I, and 




v'^(x) = STEER*^ = — -t: = , while and C2 are crowd controls. 




Thus, the crowd kinematics ([23]) expands into the matrix ODE: 
X \ I cos'' 6 \ 

y = drive'' • + STEER*" • C2 = siu'' 6 ■ 4' + I \ ■ . (24) 
9 J \ J 

A 3D simulation of random, dissipative crowd kinematics ((23|)-([24ll of 120 pcnguin-like S* £^(2)— robots, 
developed in C-|-+/DirX is presented in Figure [21 

^The above computations could instead be done in other available packages, such as Maple, by suitably translating 
the provided example code. 
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Figure 2: Driving and steering random £^(2)— dynamics of 120 penguin-like robots (with embed- 
ded collision-detection). Compare with [5]. 

The dissipative crowd kinematics ((23l) ~ ([24)) obeys the set of n-tuple integral rules of motion 
that are similar (though slightly simpler) to the above rules of the car kinematics, including the 
following derived vector-fields: 

wriggle''' = [steer''', drive'^] = [v*',u'^] and slide'' = [drive'' ,wriggle''] = [[u'', v*-'], u'']. 

Thus, controlled by the two vector controls and Cj, the crowd can form the vector-fields: 
DRIVE = u'', STEER = v'', WRIGGLE = [v'',u''], and SLIDE = [[u''', v''] , u''] , allowing it to move 
anywhere within its configuration manifold M given by (jl7[) . Solution of the dissipative crowd 
kinematics ([23|) -([24 l) defines the dissipative crowd flow, (f)^ . 

Now, the general CV-CJC crowd behavior can be defined as a amalgamate flow (behavior- 
Lagrangian flow, , plus dissipative kinematic flow, (f)^) on the crowd manifold M deflned by 

Ct=<j>i ■.t^{M{t),g{t)), 

which is a one-parameter family of honieomorphic (topologically equivalent) Riemannian mani- 
fold^ {M,g = gij), parameterized by a 'time' parameter t. That is, Ct can be used for describing 

^Proper difTcrcntiation of vector and tensor fields on a smooth Riemannian manifold (like the crowd 3n— manifold 
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smooth deformations of the crowd manifold M over time. The manifold family {M{t), g{t)) at 
time t determines the manifold family {M{t + dt),g{t + dt)) at an infinitesimal time t + dt into the 
future, according to some presecribed geometric flow, like the celebrated Ricci flow [^[^[^[21] 
(that was an instrument for a proof of a 100-year old Poincare conjecture), 



dtg,j{t) = -2R,,{t), (25) 
where Rij is the Ricci curvature tensor of the crowd manifold M and dtg[t) is defined as 

^^*.W-4(t):=lim^^^±§^. (26) 
dt dt^o dt 



3.4 Aggregate behavioral— compositional dynamics (AV) 

To formally develop the meso-level aggregate behavioral-compositional dynamics {AT>), we start 
with the crowd path integral ([201) , which can be redefined if wc Wick-rotate the time variable t to 
imaginary values, t t — it, thereby transforming the Lorentzian path integral in real time into 
the Euclidean path integral in imaginary time. Furthermore, if we rectify the time axis back to 

M) is performed using the Levi-Civita covariant derivative (see, e.g., 1411 144p . Formally, let M be a Riemannian 
A'^— manifold with the tangent bundle TM and a local coordinate system {a:*}^^ defined in an open set U C M. The 
covariant derivative operator, Vjsc : C°°{TM) —> C°°{TM), is the unique linear map such that for any vector-fields 
X, Y, Z, constant c, and scalar function / the following properties are valid: 

Vx+cy = Vx + cVy, yx{Y + fZ) = VxY + iXf)Z + f\'xZ, VxV - VyX = [X, y], 

where [X, Y] is the Lie bracket of X and Y. In local coordinates, the metric g is defined for any orthonormal basis 
{di = d/dx') in U G M by g^j = g{di,dj) = 5ij, dkgij = 0. Then the affine Levi-Civita connection is defined on 
M by 

\7g^dj = r*j9fc, where FJ^ = (pi9jl + djgu — digij) are the Christoffel symbols. 

Now, using the covariant derivative operator Vx we can define the Riemann curvature (3, 1)— tensor *Rm by 

^Rm(X, Y)Z = Vx^yZ - Vy^xZ - Vyx,Y]Z, 

which measures the curvature of the manifold by expressing how noncommutative covariant differentiation is. The 
(3, 1)— components R\jf. of 5Hm are defined inU d M by 

^ra{d,,dj)ak = R\,^du or iJ^.fc = a^rj^ - 5,rjfe + Tf, F^ - T^r^™- 

Also, the Riemann (4,0)— tensor Rijkl = dlmR^i^ is defined as the g— based inner product on M, 

Rijki = (iHm(9,,9j)9fe,0;) . 
The first and second Bianchi identities for the Riemann (4,0)— tensor Rijkl hold, 

Rijkl + Rjkil + Rkijl = 0, ^iRjklm + ^ j Rkilm + ^ kRijlm = 0, 

while the twice contracted second Bianchi identity reads; 2\/jRij = V^i?. 
The (0,2) Ricci tensor IHc is the trace of the Riemann (3, 1)— tensor *Hm, 

5Kc(y, Z) + tr(X ^ lKm(X,y)Z), so that 9i<:(X,Y) = gCCfim{di, X)di,Y), 

Its components Rjk = ^Hc {dj, dk) are given in (/ C M by the contraction 

Rjk = R\jk^ °r Rjk = diT)k - dk'^)i + ^Ini^Tk - r^fcr™. 

Finally, the scalar curvature R is the trace of the Ricci tensor JHc, given inU <Z M by: R = g^^ Rij. 
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the real line, we get the adaptive SFT~partition function as our proposed ^27-model: 

(Physical Action | Mental Preparation) ad = / T)[w, x, y] e-^[^'f'*-*il. (27) 

The adaptive ^2?-transition amplitude (Physical Action | Mental Preparation) ad as defined 
by the SFT-partition function (j27[) is a general model of aMarkov stochastic process. Recall 
that Markov process is a random process characterized by a lack of memory, i.e., the statistical 
properties of the immediate future are uniquely determined by the present, regardless of the past 
(see, e.g. [TSl SI])- The iV— dimensional Markov process can be defined by the Ito stochastic 
diS'erential equation, 

dxiit) = A,[x\t),t]dt + B,^[x'{t),t]dW^{t), (28) 
x'iO) = x,o, (z,j = l,...,iV) (29) 

or corresponding Ito stochastic integral equation 

x'{t)=x'{0)+ [ dsA^[x'{s),s]+ [ dW^s)B,j[x\s),s], (30) 



in which x^{t) is the variable of interest, the vector Ai[x{t), t] denotes deterministic drift, the matrix 
Bij[x{t),t\ represents continuous stochastic diffusion fluctuations, and W^{t) is an N— variable 
Wiener process (i.e., generalized Brownian motion [15j ) and 

dW^{t) = W^{t + dt) - W^{t). 

The two Ito equations (|29p - (|30p are equivalent to the general Chapman-Kolmogorov probability 
equation (see equation pip below). There are three well known special cases of the Chapman- 
Kolmogorov equation (see [T^ ) : 

1. When both Bij[x{t),t] and W{t) are zero, i.e., in the case of pure deterministic motion, it 
reduces to the Liouville equation 

dtPix', t'\x\ t") = - E ^ {ArW). t] P{x\t'\x", t")} . 

i 

2. When only W{t) is zero, it reduces to the Fokker-Planck equation 

d 

dx^ 

+ lT.7^iB.^i-it)^t]p{x',tv'X)}. 



dtP{x',t'\x",t") = -Y^l-{AMt),t]P{x',t'\x",t")} 



2 dx^dxi 

3. When both Ai[x{t),t] and Bij[x{t),t\ are zero, i.e., the state-space consists of integers only, 
it reduces to the Master equation of discontinuous jumps 



dtPix',t'\x",t") = / dxW{x'\x",t)P{x',t'\x",t")- / dxWix"\x',t)P{x',t'\x",t") 
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The Markov assumption can now be formulated in terms of the conditional probabilities 
P{x^,ti): if the times ti increase from right to left, the conditional probability is determined 
entirely by the knowledge of the most recent condition. Markov process is generated by a set of 
conditional probabilities whose probability-density P = P{x' , t'\x" , t") evolution obeys the general 
Chapman-Kolmogorov integro-differential equation 

i ij 

including deterministic drift, diffusion fluctuations and discontinuous jumps (given respectively 
in the first, second and third terms on the r.h.s.). This general Chapman-Kolmogorov integro- 
differential equation (|3ip. with its conditional probability density evolution, P = P(x', t"), 
is represented by our SFT-partition function (|27p . 

Furthermore, discretization of the adaptive SFT-partition function (|27p gives the standard 
partition function 

Z = ^e-"'^^VT^ (31) 

3 

where is the motion energy eigenvalue (reflecting each possible motivational energetic state), 
T is the temperature-like environmental control parameter, and the sum runs over all ID energy 
eigenstates (labelled by the index j). From ((3T|) . we can calculate the transition entropy, as 
S* = fcs InZ (see the next section). 

4 Entropy, chaos and phase transitions in the crowd 
manifold 

Recall that nonequilibrium phase transitions [T71 [T51 [121 HOI HI] are phenomena which bring about 
qualitative physical changes at the macroscopic level in presence of the same microscopic forces 
acting among the constituents of a system. In this section we extend the CD formalism to incorpo- 
rate both algorithmic and geometrical entropy as well as dynamical chaos [5H [5H [5D] between 
the entropy-growing phase of Mental Preparation and the entropy-conserving phase of Physical 
Action, together with the associated topological phase transitions. 

4.1 Algorithmic entropy 

The Boltzmann and Shannon (hence also Gibbs entropy, which is Shannon entropy scaled by k In 2, 
where k is the Bolzmann constant) entropy definitions involve the notion of ensembles. Membership 
of microscopic states in ensembles defines the probability density function that underpins the 
entropy function; the result is that the entropy of a definite and completely known microscopic 
state is precisely zero. Bolzmann entropy defines the probabilistic model of the system by effectively 
discarding part of the information about the system, while the Shannon entropy is concerned with 
measuring the ignorance of the observer - the amount of missing information - about the system. 
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Zurck proposed a new physical entropy measure that can be apphed to individual microscopic 
system states and does not use the ensemble structure. This is based on the notion of a fixed 
individually random object provided by Algorithmic Information Theory and Kolmogorov Com- 
plexity: put simply, the randomness K{x) of a binary string x is the length in terms of number of 
bits of the smallest program p on a universal computer that can produce x. 

While this is the basic idea, there are some important technical details involved with this defini- 
tion. The randomness definition uses the prefix complexity K(.) rather than the older Kolmogorov 
complexity measure C(.): the prefix complexity K{x\y) of x given y is the Kolmogorov complexity 
C^^{x\y) = min{p | x = (pu{{y,p))} (with the convention that C^^{x\y) = c» if there is no suchp) 
that is taken with respect to a reference universal partial recursive function that is a universal 
prefix function. Then the prefix complexity K{x) of x is just K{x\e) where e is the empty string. 
A partial recursive prefix function ^ : M — > N is a partial recursive function such that if ^{p) < oo 
and (f>{q) < oo then p is not a proper prefix of q: that is, we restrict the complexity definition to 
a set of strings (which are descriptions of effective procedures) such that none is a proper prefix 
of any other. In this way, all effective procedure descriptions are self-delimiting: the total length 
of the description is given within the description itself. A universal prefix function (/)„is a prefix 
function such that \/n G N <pu{{y, {n,p))) = <pn{{y,p))j where <pn is numbered n according to some 
Godel numbering of the partial recursive functions; that is, a universal prefix function is a partial 
reciirsivc function that simulates any partial recursive function. Here, (x, y) stands for a total 
recusive one-one mapping from NxN into N, {xi,X2, ■ ■ ■ ,Xn) = {xi, {x2, ■ ■ ■ ,Xn)), N is the set of 
natural numbers, and M = {0, l}*is the set of all binary strings. 

This notion of entropy circumvents the use of probability to give a concept of entropy that can be 
applied to a fully specified macroscopic state: the algorithmic randomness of the state is the length 
of the shortest possible effective description of it. To illustrate, suppose for the moment that the 
set of microscopic states is countably infinite, with each state identified with some natural number. 
It is known that the discrete version of the Gibbs entropy (and hence of Shannon's entropy) and 
the algorithmic entropy are asymptotically consistent under mild assumptions. Consider a system 
with a countably infinite set of microscopic states X supporting a probability density function P{.) 
so that P{x) is the probability that the system is in microscopic state x £ X. Then the Gibbs 

entropy is Sg{P) = — (fcln2) ^ P(a;) logP(a;) (which is Shannon's information-theoretic entropy 

xex 

H{P) scaled by A; In 2). Supposing that P(.) is recursive, then Sg{P) = (A; In 2) P{x)K{x)+C, 

xex 

where is a constant depending only on the choice of the reference universal prefix function ^. 
Hence, as a mcasm-e of entropy, the function K{.) manifests the same kind of behavior as Shannon's 
and Gibbs entropy measures. 

Zurek's proposal was of a new physical entropy measure that includes contributions from both 
the randomness of a state and ignorance about it. Assume now that we have determined the 
macroscopic parameters of the system, and encode this as a string - which can always be converted 
into an equivalent binary string, which is just a natural number under a standard encoding. It is 
standard to denote the binary string and its corresponding natural mmiber interchangeably; here 
let X be the encoded macroscopic parameters. Zurek's definition of algorithmic entropy of the 
macroscopic state is then K{x) + H^, where = 5_B(a;)/(fcln2), where Sb{x) is the Bolzmann 
entropy of the system constrained by x and k is Bolzmann's constant; the physical version of the 
algorithmic entropy is therefore defined as Sa{x) = {klii2){K{x) + Hx). Here represents the 
level of ignorance about the microscopic state, given the parameter set x; it can decrease towards 
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zero as knowledge about the state of the system increases, at which point the algorithmic entropy 
reduces to the Bolzmann entropy. 

4.2 Ricci flow and Perelman entropy— action on the crowd manifold 

Recall that the inertial metric crowd flow, : i i— > (M(t), g(t)) on the crowd 3n— manifold (fT7| is 
a one-parameter family of homeomorphic Riemannian manifolds (M, 17), evolving by the Ricci flow 

Now, given a smooth scalar function u : Af — > R on the Riemannian crowd 3n— manifold M, 
its Laplacian operator A is locally deflned as 

Au = g'^V.VjU, 

where Vi is the covariant derivative (or, Levi-Civita connection). We say that a smooth function 
u : M X [0, T) E, where T e (0, cxo], is a solution to the heat equation on M if 

dtu = Au. (32) 

One of the most important properties satisfied by the heat equation is the maximum principle, 
which says that for any smooth solution to the heat equation, whatever point-wise bounds hold at 
t = also hold for t > [Hj . This property exhibits the smoothing behavior of the heat diffusion 
(021) on M. 

Closely related to the heat diffusion ([5^ is the (the Fields medal winning) Perelman entropy- 
action functional, which is on a 3n— manifold M with a Riemannian metric gij and a (temperature- 
like) scalar function / given by [64j 

£ = f (R+\Vf\^)e-fdn (33) 

JM 

where R is the scalar Riemann curvature on M, while dfi is the volume 3?!— form on M, defined as 



dfj. = y det(g,y) dx^ A dx"^ A ... A dx^". (34) 

During the Ricci fiow ((25)) " ([26l) on the crowd manifold (|17p . that is, during the inertial metric 
crowd flow, Ct'.t^ {M{t), g{t)), the Perelman entropy functional ([55)1 evolves as 

dt£ = 2 J \R,j + V,;Vj/pe"^ d^l. (35) 

Now, the crowd breathers are solitonic crowd behaviors, which could be given by localized 
periodic solutions of some nonlinear soliton PDEs, including the exactly solvable sine-Gordon 
equation and the focusing nonlinear Schrodinger equation. In particular, the time-dependent crowd 
inertial metric gij{t), evolving by the Ricci flow g{t) given by (P5)) - (PS1) on the crowd 3ri— manifold 
M is the Ricci crowd breather, if for some ti < t2 and a > the metrics agij{ti) and gij{t2) 
differ only by a diffeomorphism; the cases a — \,a<l,a> \ correspond to steady, shrinking and 
expanding crowd breathers, respectively. Trivial crowd breathers, for which the metrics gij{ti) 
and gij{t2) on A'l differ only by diffeomorphism and scaling for each pair of ti and t2, are the 
crowd Ricci solitons. Thus, if we consider the Ricci flow (|25p - (|26p as a biodynamical system on 
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the space of Riemannian metrics modulo difFeomorphism and scaling, then crowd breathers and 
solitons correspond to periodic orbits and fixed points respectively. At each time the Ricci soliton 
metric satisfies on M an equation of the form [51] 

Rij +cg,j +Vibj +Vjbi = 0, 

where c is a number and bi is a 1-form; in particular, when bi — ^V^a for some function a on M, 
we get a gradient Ricci soliton. 

Define X{gij) ~ ini £{gij , f), where infimum is taken over all smooth /, satisfying 

/ e-^dfi = 1. (36) 
Jm 

\{gij) is the lowest eigenvalue of the operator — 4A + R. Then the entropy evolution formula pS]) 
imphes that \{gij{t)) is non-decreasing in i, and moreover, if \{ti) = A(t2), then for t G [ii,t2] we 
have Rij + ViVjf = for / which minimizes £" on M [51]. Therefore, a steady breather on M is 
necessarily a steady soliton. 

If we define the conjugate heat operator on M as 

□* = -d/dt -A + R 

then we have the conjugate heat equation: \I\*u = 0. 

The entropy functional ([551) i'' nondecreasing under the coupled Ricci-diffusion flow on M [IH] 

dtg^j = -2R,j, dtu = -Au +^u~ (37) 

where the second equation ensures J^^ v?d^ = 1, to be preserved by the Ricci flow g(t) on M . If 

we deflne u = , then (|37p is equivalent to /—evolution equation on M (the nonlinear backward 
heat equation), 

ft/ = -A/+|V/|2-i?, 

which instead preserves ([55)1 . The coupled Ricci-diffusion flow (1571) is the most general biodynamic 
model of the crowd reaction-diffusion processes on M . In a recent study [1] this general model 
has been implemented for modelling a generic perception-action cycle with applications to robot 
navigation in the form of a dynamical grid. 

Perelman's functional £ is analogous to negative thermodynamic entropy [64 . Recall that 
thermodynamic partition function for a generic canonical ensemble at temperature is given by 

j e-^^diuiE), (38) 

where uj{E) is a 'density measure', which does not depend on f3. From it, the average energy is given 
by (E) = -dp In Z, the entropy is S = (3 {E) + In Z, and the fluctuation is a = {{E - {E)Y) = 
d/32 In Z. 

If we now fix a closed 3rt— manifold M with a probability measure m and a metric gijir) that 
depends on the temperature r, then according to equation 

^rg^J = 2(i?„ + V,Vj/), 
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the partition function (|38p is given by 



lnZ = J {-/ + -) dm. (39) 

From ((391) we get (see [64]) 

Tl f 

— )dm, S^- {t{R + I V/p) + f-n) dm, 

— (yfy P dm, where dm = udV, u = (Attt)^^ e^-^ . 

From the above formulas, we see that the fluctuation a is nonnegative; it vanishes only on a 
gradient shrinking soliton. (E) is nonnegative as well, whenever the flow exists for all sufflciently 
small T > 0. Furthermore, if the heat function u: (a) tends to a (5— function as r ^ 0, or (b) 
is a limit of a sequence of partial heat functions Ui, such that each Ui tends to a (5— function as 
T Ti > 0, and Ti — > 0, then the entropy S is also nonnegative. In case (a), all the quantities 
(E) ,S,a tend to zero as r — > 0, while in case (b), which may be interesting if gijir) becomes 
singular at r = 0, the entropy S may tend to a positive limit. 



{E) = -t' I (i?+|V/p- 

JM 

a = 2T^ [ + V.Vj/ - 

JM 



4.3 Chaotic inter-phase in crowd dynamics induced by its Riemannian 
geometry change 

Recall that CD transition map ([5]) is defined by the chaotic crowd phase-transition amplitude 

/ dtS=Q dtS>0 \ f , 

/PHYS. ACTION CHAOS MENTAL PREP.\ := J V[x]e'^^''\ 

where we expect the inter-phase chaotic behavior (see [45l ) . To show that this chaotic inter-phase is 
caused by the change in Riemannian geometry of the crowd 3n— manifold Af , we will first simplify 
the CV action functional (jl8|) as 

A[x] ^11"" [9^J i'i' - Vix, x)] dt, (40) 
with the associated standard Hamiltonian, corresponding to the amalgamate version of (|14p . 

N 

2' 



" 1 

^b'2:) =^^P? + ^(a:,i), (41) 



where Pi are the SE(2)-momenta, canonically conjugate to the individual agents' SE(2)-coordinates 
x% (i = 1, 3n). Biodynamics of systems with action ([^01) and Hamiltonian (HT1) are given by the 
set of geodesic equations [HJ Si] 

^ + r;..^^ = o, (42) 

ds^ ds ds 

where F*^ are the Christoffel symbols of the affine Levi-Civita connection of the Riemannian CD 
manifold M . In this geometrical framework, the instability of the trajectories is the instability of 
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the geodesies, and it is completely determined by the curvature properties of the CV manifold M 
according to the Jacobi equation of geodesic deviation [411 Hll 



^+^.7c™^J ^ = 0, (43) 

whose solution J, usually called Jacobi variation field, locally measures the distance between nearby 
geodesies; D/ds stands for the covariant derivative along a geodesic and j]^^ are the components 
of the Riemann curvature tensor of the CD manifold M . 

The relevant part of the Jacobi equation (|43)) is given by the tangent dynamics equation [8l [5] 

J^+i?Vo-/' = 0, (*,A: = l,...,3n), (44) 

where the only non-vanishing components of the curvature tensor of the CD manifold M are 

= d^V/dx^dxK (45) 

The tangent dynamics equation (|44p can be used to define Lyapunov exponents in dynamical 
systems given by the Riemannian action (|40p and Hamiltonian ()4ip . using the formula [7] 

Ai = hm l/2nog(Ad[jf(t) + J^{t)yMl^[J^{Q) + 7^(0)]). (46) 

t — >oo 

Lyapunov exponents measure the strength of dynamical chaos in the crowd behavior dynamics. 
The sum of positive Lyapunov exponents defines the Kolmogorov-Sinai entropy. 



4.4 Crowd nonequilibrium phase transitions induced by manifold 
topology change 

Now, to relate these results to topological phase transitions within the CD manifold M given by 
p7|) . recall that any two high-dimensional manifolds My and M^' have the same topology if they 
can be continuously and differentiably deformed into one another, that is if they are diffeomorphic. 
Thus by topology change the 'loss of diffeomorphicity' is meant [67]. In this respect, the so- 
called topological theorem 43] says that non-analyticity is the 'shadow' of a more fundamental 
phenomenon occurring in the system's configuration manifold (in our case the CD manifold): a 
topology change within the family of equipotcntial hypersurfaces 

M„ = {{x\...,x^") gR3"| y(x\...,x=^") 

where V and are the microscopic interaction potential and coordinates respectively. This topo- 
logical approach to PTs stems from the numerical study of the dynamical counterpart of phase 
transitions, and precisely from the observation of discontinuous or cuspy patterns displayed by the 
largest Lyapunov exponent Ai at the transition energy |7]. Lyapunov exponents cannot be mea- 
sured in laboratory experiments, at variance with thermodynamic observables, thus, being genuine 
dynamical observables they are only be estimated in numerical simulations of the microscopic dy- 
namics. If there are critical points of V in configuration space, that is points Xc — [xi, . . . ,X3„] 
such that Vy(x)|^^^ = 0, according to the Morse Lemma [32], in the neighborhood of any critical 
point Xc there always exists a coordinate system x{t) — [x^{t), ...,x^'^{t)] for which |7] 

V{x) = Vixc) ~xl xl+ + • • • + xl„, (47) 
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where k is the index of the critical point, i.e., the number of negative eigenvalues of the Hessian of 
the potential energy V. In the neighborhood of a critical point of the CP-manifold M, equation 
([TT)) yields the simplified form of d^V/dx^dx^ = ±5y, giving j unstable directions that 

contribute to the exponential growth of the norm of the tangent vector J. 

This means that the strength of dynamical chaos within the C2?-manifold M, measured by 
the largest Lyapunov exponent Ai given by is affected by the existence of critical points Xc 
of the potential energy V{x). However, as V(x) is bounded below, it is a good Morse function, 
with no vanishing eigenvalues of its Hessian matrix. According to Morse theory [32] , the existence 
of critical points of V is associated with topology changes of the hypersurfaces {M„}„gR. The 
topology change of the {Mi,}„gH at some Vc is a necessary condition for a phase transition to take 
place at the corresponding energy value [13] . The topology changes implied here are those described 
within the framework of Morse theory through 'attachment of handles' [32] to the CD-manifold 
M. 

In our path-integral language this means that suitable topology changes of equipotential sub- 
manifolds of the CP-manifold M can entail thermodynamic-like phase transitions [TT] UHl [19], 
according to the general formula: 



(phase out | phase in) := / 2?[w$] e'^I*' 

<y top— ch 



The statistical behavior of the crowd biodynamics system with the action functional (HDl) and the 
Hamiltonian (|4ip is encompassed, in the canonical ensemble, by its partition function, given by 
the Hamiltonian path integral [44] 

Z3n= f V[p]V[x]exp{i f [p,x' ~H{p,x)]dr}, (48) 



dxlT)dp{T) 
2^ ■ 



J top — ch 

where we have used the shorthand notation 



/ v[p]v[x] ^ / n 

•^toD— ch -/ „ 



/ top— ch 

The path integral (UHl) be calculated as the partition function 



/oil / \ ^ r 

n dp, dx'e-^"^^'-^ = i^^j J n dx'^ 

^^dve-^^ I J^, (49) 



Jo Jm^ W^v 

where the last term is written using the so-called co-area formula [lOj . and v labels the equipo- 
tential hypersurfaces M„ of the CD manifold M, 

M„ = {{x\ 2)3") e R^''\V{x\ . . . ,x3") = v}. 

Equation (jl^ shows that the relevant statistical information is contained in the canonical config- 
urational partition function 

Zg- jY[dx^Vix)e-^^^^\ 
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Note that Z^^^ is decomposed, in the last term of (|^^ . into an infinite summation of geometric 
integrals, 

/ da/\\VVl 

defined on the {M„}^gR. Once the microscopic interaction potential V{x) is given, the configura- 
tion space of the system is automatically foliated into the family {M^}y^s. of these equipotential 
hypersurfaces. Now, from standard statistical mechanical arguments we know that, at any given 
value of the inverse temperature /3, the larger the number 3n, the closer to My = are the 
microstates that significantly contribute to the averages, computed through Z^niP), of thermody- 
namic observables. The hypersurface is the one associated with 

u, = {Z^J-'ll[dx^V{x)e-f'^^^\ 

the average potential energy computed at a given /3. Thus, at any f3, if 3n is very large the 
effective support of the canonical measure shrinks very close to a single My = . Hence, the 
basic origin of a phase transition lies in a suitable topology change of the {My}, occurring at some 
Vc [12] • This topology change induces the singular behavior of the thermodynamic observables at 
a phase transition, ft is conjectured that the counterpart of a phase transition is a breaking of 
diffeomorphicity among the surfaces My , it is appropriate to choose a diffeomorphism invariant to 
probe if and how the topology of the My changes as a function of v. Fortunately, such a topological 
invariant exists, the Euler characteristic of the crowd manifold M, defined by [HI |44] 

3n 

x(M)=^(-l)'^6fe(M), (50) 

where the Betti numbers bk{M) are diffeomorphism invariants are the dimensions of the de 
Rham's cohomology groups 7J'=(M;M); therefore the bk are integers). This homological formula 
can be simplified by the use of the Gauss-Bonnet theorem, that relates xi^) with the total 
Gauss-Kronecker curvature Kg of the C2?-manifold M given by [HI |SD| 

x(M) = / Kcdfi, where dfi is given by ([M]) . 
Jm 

5 Conclusion 

Our understanding of crowd dynamics is presently limited in important ways; in particular, the 
lack of a geometrically predictive theory of crowd behavior restricts the ability for authorities to 
intervene appropriately, or even to recognize when such intervention is needed. This is not merely 
an idle theoretical investigation: given increasing population sizes and thus increasing opportunity 
for the formation of large congregations of people, death and injury due to trampling and crushing 
- even within crowds that have not formed under common malicious intent - is a growing concern 
among police, military and emergency services. This paper represents a contribution towards the 
understanding of crowd behavior for the purpose of better informing decision-makers about the 
dangers and likely consequences of different intervention strategies in particular circumstances. 
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In this paper, we have proposed an entropic geometrical model of crowd dynamics, with dis- 
sipative kinematics, that operates across macro-, micro- and meso-levels. This proposition is 
motivated by the need to explain the dynamics of crowds across these levels simultaneously. We 
contend that only by doing this can we expect to adequately characterize the geometrical properties 
of crowds with respect to regimes of behavior and the changes of state that mark the boundaries 
between such regimes. 

In pursuing this idea, we have set aside traditional assumptions with respect to the separation 
of mind and body. Furthermore, we have attempted to transcend the long-running debate between 
contagion and convergence theories of crowd behavior with our multi-layered approach: rather than 
representing a reduction of the whole into parts or the emergence of the whole from the parts, our 
approach is build on the supposition that the direction of logical implication can and does flow 
in both directions simultaneously. We refer to this third alternative, which effectively unifies the 
other two, as behavioral composition. 

The most natural statistical descriptor is crowd entropy, which satisfies the extended second 
thermodynamics law applicable to open systems comprised of many components. Similarities 
between the configuration manifolds of individual (micro-level) and crowds (macro-level) motivate 
our claim that goal-directed movement operates under entropy conservation, while natural crowd 
dynamics operates under monotonically increasing entropy functions. Of particular interest is 
what happens between these distinct topological phases: the phase transition is marked by chaotic 
movement. 

We contend that this approach provides a basis on which one can build a geometrically predic- 
tive model of crowd behavior dynamics - over and above the existing approaches, which are largely 
explanatory. The current paper develops an entropy formulation of crowd dynamics as a three-step 
process involving individual and collective behavior dynamics, and - crucially - non-equilibrium 
phase transitions whereby the forces operating at the microscopic level result in geometrical change 
at the macroscopic level. We have incorporated both geometrical and algorithmic notions of en- 
tropy as well as chaos in studying the topological phase transition between the entropy conservation 
of physical action and the entropy increase during internal, action preparation stages of behavior. 
Given these formulations, future research can focus on: (i) crowd simulations in 3D graphics envi- 
ronments, (ii) motivated cognition underpinning crowd dynamics and (iii) mechanisms of abrupt 
change in crowd behavior, such as crowd turbulence and flashpoints. 
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6 Appendix 

General nonlinear stochastic dynamics, developed in a framework of Feynman path integrals, have 
recently [46] been applied to Lewinian field-theoretic psychodynamics [58|, resulting in the devel- 
opment of a new concept of Life-Space Foam (LSF) as a natural medium for motivational and 
cognitive psychodynamics. According to the LSF-formalism, the classic Lewinian life space can be 
macroscopically represented as a smooth manifold with steady force-fields and behavioral paths. 
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while at the microscopic level it is more realistically represented as a collection of wildly fluctuating 
force-fields, (loco)motion paths and local geometries (and topologies with holes). 

A set of least-action principles is used to model the smoothness of global, macro-level LSF 
paths, fields and geometry, according to the following prescription. The action S'[$], with dimen- 
sions of Energy x Time = Effort and depending on macroscopic paths, fields and geometries 
(commonly denoted by an abstract field symbol (f>*) is defined as a temporal integral from the 
initial time instant ti„i to the final time instant 



sm = m dt, (51) 

with Lagrangian density given by 

£[$] = J d"x/:($„a,.$^), 

where the integral is taken over all n coordinates = x^{t) of the LSF, and d^i^^ are time and 
space partial derivatives of the variables over coordinates. The standard least action principle 

= 0, (52) 

gives, in the form of the so-called Euler-Lagrangian equations, a shortest (loco)motion path, an 
extreme force-field, and a life-space geometry of minimal curvature (and without holes). In this 
way, we have obtained macro-objects in the global LSF: a single path described by Newtonian-like 
equation of motion, a single force-field described by Maxwellian-like field equations, and a single 
obstacle-free Riemannian geometry (with global topology without holes). 

To model the corresponding local, micro-level LSF structures of rapidly fluctuating MD & CD, 
an adaptive path integral is formulated, defining a multi-phase and multi-path (multi-field and 
multi-geometry) transition amplitude from the motivational state of Intention to the cognitive 
state of Action, 

{Action\Intention) total ■■= ^ X»[w$] e''^'*' , (53) 

where the Lebesgue integration is performed over all continuous $con ~ paths+fields+georiietries, 
while summation is performed over all discrete processes and regional topologies <1>^^^ . The symbolic 
differential I?[w<I>] in the general path integral represents an adaptive path measure, defined 
as a weighted product 



N 

I?[w<i>] = lim TTiy^d^'g, {i = 1, ...,n = con + dis). (54) 

TV ^ ^ 



The adaptive path integral (|53p -(|54 p represents an oo— dimensional neural network, with weights 
w updating by the general rule [49j 

new value(t + 1) = old value{t) + innovation{t) . 
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